Conditional probabilities with Dirac observables and the problem of time in quantum 

gravity 
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We combine the "evolving constants" approach to the construction of observables in canonical 
quantum gravity with the Page-Wootters formulation of quantum mechanics with a relational time 
for generally covariant systems. This overcomes the objections levied by Kuchaf against the latter 
formalism. The construction is formulated entirely in terms of Dirac observables, avoiding in all 
cases the physical observation of quantities that do not belong in the physical Hilbert space. We work 
out explicitly the example of the parameterized particle, including the calculation of the propagator. 
The resulting theory also predicts a fundamental mechanism of decoherence. 



In generally covariant systems, like general relativity, 
when one works out the canonical formulation the Hamil- 
tonian is a constraint, i.e. it vanishes identically. This 
implies that the parameter that usually plays the role of 
time in canonical formulations is not adequate to describe 
the dynamics of the system. This constitutes one of the 
aspects of the "problem of time" for generally covariant 
systems (see Kuchaf [lj for a good review). Page and 
Wootters proposed an approach to deal with this issue 
0. The proposal consists in building a quantum theory 
of the system of interest by promoting all variables of the 
system to quantum operators and then choosing one of 
the variables to be a "clock" and computing conditional 
probabilities for the other variables to take certain values 
when the "clock" takes a given value. This proposal ran 
into technical difficulties when applied in detail to con- 
strained systems, as emphasized by Kuchaf [l[ . Basically, 
the problem consists in what to choose as the variables 
to be promoted to operators, in particular which one to 
choose as a "clock" . In a constrained system the physi- 
cally observable variables are those that have vanishing 
Poisson brackets with the constraints (this implies they 
are invariant under the symmetries of the theory, they 
are "gauge invariant" ) . However, if one of the constraints 
is the Hamiltonian, then quantities that have vanishing 
Poisson brackets with it do not evolve and therefore are 
poor candidates for being clocks. Page and Wootters 
tried to circumvent this by considering "kinematical" 
variables (that do not have vanishing Poisson brackets 
with the constraints and therefore appear "to evolve"). 
But this brings about other problems. Such variables can 
be promoted to quantum operators acting on the space of 
wavefunctions that are not necessarily annihilated by the 
constraints ("kinematical Hilbert space"). Within such 
space, the states that are annihilated by the constraints 
are usually distributional (at least in simple examples 
that can be worked out explicitly). Therefore they may 
not admit a probabilistic interpretation. Kuchaf showed 
by analyzing the example of a parameterized particle, 
that these issues had as a consequence the construction 



of propagators that "do not propagate" and therefore the 
resulting quantum theory is not realistic. 

Here we would like to revisit the Page-Wootters con- 
struction but using a different set of physical quantities. 
The quantities we will choose are relational Dirac observ- 
ables such as the "evolving constants" introduced in Q 
(an idea that goes back to DeWitt, Bergmann and Ein- 
stein). The proposal can be summarized as follows. In 
a totally constrained theory, the values of fields are not 
physically observable. On the other hand, if one chooses a 
one-parameter family of observables such that their value 
coincides with the value of a dynamical variable when the 
parameter takes the value of another dynamical variable, 
which one uses to characterize the evolution, such ob- 
servables can be used in the Page-Wootters construction. 
They have the advantage that there is a sense in which 
they "evolve" . That is, unlike the proposal of Rovelli, 
we will not consider the "parameter" to be the physical 
time, but we will use it to make sense of the conditional 
probabilities that arise in the Page-Wootters formulation 
when one introduces a real quantum clock. In fact at the 
end of the day the parameter drops out from the formu- 
lation, and one integrates over all possible values of it 
(if one has more than one constraint one needs to intro- 
duce more parameters). Therefore, one does not need 
to observe any dynamical variable that is not quantum 
mechanical or is not a Dirac observable. We will show in 
an example that this construction can be carried out in 
detail. 

Let us start by defining the evolving constants in a 
classical theory. Following Rovelli we consider a totally 
constrained system with a phase space q l ,p l - We now 
pick a parameter we call t. We are interested in defining 
a one-parameter family of Dirac observables that repro- 
duces the value of one of the dynamical variables, for 
instance q l , when another variable playing the role of a 
clock takes the value t. For concreteness, we can choose 
q 1 to play the role of a clock. We denote the one param- 
eter family of Dirac observables Q l (t) — Q l (t,q n ,p n ). 
These have vanishing Poisson bracket with the Hamilto- 
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nian constraint, {Q l {t, q n : p n ) : C(q n ,p n )} = 0. They are 
also such that if one evaluates Q l (t, t,q 2 . . . q n ,p n ) = q l . 
(We refer to the observables as Q l for simplicity, they can 
include momenta as well, but they must have vanishing 
Poisson bracket with the clock variable, an assumption 
that may be relaxed with further elaboration) . 

We then proceed to quantize the variables of the prob- 
lem. Namely, we promote all Q*(t)'s and P l (t)'s to 
quantum self-adjoint operators acting on a Hilbert space 
of wavefunction ^(q l ) that are annihilated by the con- 
straints. The variable t will remain classical. In realistic 
situations, like general relativity, this is convenient since 
usual choices of "time" are given by global variables that 
are not easily associated with a quantum operator to be- 
gin with. The restriction to self adjoint operators lim- 
its importantly the choices of possible parameters t, as 
was discussed in In particular if one docs not insist 
on self-adjointness one runs into problems in the defi- 
nition of conditional probabilities (related to the "false 
tracks" discussed in [5j). If the quantization can be ac- 
complished, then one can introduce a basis of eigenstates 
(parameterized by t) of the evolving constants and intro- 
duce projectors that materialize the physical properties 
associated to each of the evolving constants. 

Having quantized the evolving constants, we will 
choose one of the variables Q l (t) to be a quantum clock, 
and we will call it T{t). We then introduce the condi- 
tional probability as in j(| @] 



P(Q'=Q' \t=t ) = lim 



f_ T dm[p Q}s {t)p To {t) P p To {t) 
f_ T dm(p To (t)p) 



(i) 

where PQ' o (t) is the projector on the eigenspace associ- 
ated with the eigenvalue Q l at time t and similarly for 
PT {t). These conditional probabilities are positive and 
add to one. 

By construction the conditional probability is a gauge 
invariant quantity since the density matrix, p in the 
above expression, is assumed to be annihilated by the 
constraints, e.g. Cp = 0. Note that we are treating 
the variable t as an unobservable quantity and summing 
over all possible values of it. This picture is much more 
satisfactory than the one that emerges from considering 
evolving constants alone without the conditional proba- 
bility interpretation, since in that picture one assumes 
that a quantum variable like q , takes a definite classi- 
cal value. This would not usually be the case since q 1 
has a non-vanishing Poisson bracket with the constraint 
and on the constraint surface we expect q 1 to have infi- 
nite uncertainty. Returning to the above expression, it 
should be noted that the improper limits of integration 
may cause problems at the time of computing the con- 
ditional probabilities. This can be controlled by simply 
considering integrals in temporal domains that are much 
larger than the region of physical interest. 



It is worthwhile expanding on the meaning of the prob- 
abilities {]]) since there has been some confusion in the 
literature [8(. One may interpret that the numerator of 
([1]) is the sum of joint probabilities of O and T for all 
values of t. This would be incorrect since the events in 
different t's are not mutually exclusive. The probabil- 
ity {T]) corresponds to a physically measurable quantity, 
and that such quantity is actually the only thing one can 
expect to measure in systems where one does not have 
direct access to the "ideal" time t. The experimental 
setup we have in mind is to consider an ensemble of non- 
interacting systems with two quantum variables each to 
be measured, O and T. Each system is equipped with a 
recording device that takes a single snapshot of O and T 
at a random unknown value of the "ideal" time t. One 
takes a large number of such systems, launches them all 
in the same quantum state, "waits for a long time" , and 
concludes the experiment. The recordings taken by the 
devices are then collected and analyzed all together. One 
computes how many times n(Tj, Oj) each reading with a 
given value T = Tj,0 = Oj occurs (to simplify things, 
for the moment let us assume T, O have discrete spectra; 
for continuous spectra one would have to consider values 
in a small finite interval of the value of interest). If one 
takes each of those values n(Tj,Oj) and divides them by 
the number of systems in the ensemble, one obtains, in 
the limit of infinite systems, a joint probability P(Oj,Tj) 
that is represented by the above expression. 

We can then write the conditional probabilities that 
yield the correlation functions (propagators), namely the 
probability that having the system been observed at Q\ 
at time T\ and it will be observed at Q\ at time T 2 , 

P(Qi\T 2 ,Q\,T u p)= (2) 
f_ T dt j: r dt'Tv (P Q% ^ (t)P Q{<Tl (t')pP QlTl (O) 



lim 



f_ r dt f_ T dt'Tr (P T2 (t)P Q , >Tl (t') P P Qi , Tl (<')) 



Implicit is the use of a reduction postulate after the mea- 
surement of Q\,Ti (see 0). Notice that in particular no 
assumption about the relative ordering of the unobserv- 
able variables t and t' is needed. We will show that it 
yields the correct propagator in an example. 

Up to now recovering the correct propagator has been 
problematic in the conditional probability approach. For 
instance, Kuchaf [l| computed a similar expression us- 
ing the original Page and Wootters 0] prescription (in 
that case, however there is no t or t' and no integrals 
over them) where the quantities Q l were kinematical 
operators that did not commute with the constraint 
and showed that one obtained an incorrect propaga- 
tor. Essentially, the system did not move, the propa- 
gator being proportional to a Dirac delta function, e.g. 
<KQ| - Q{)S(T 2 - Ti). Page [9] has responded to this 
criticism by claiming that in the conditional probability 
framework one cannot compute two time probabilities. 
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We believe that the framework can indeed accommodate 
such probabilities and therefore becomes more powerful 
when formulated in terms of evolving constants and in- 
deed yields the correct propagators. 

The example we will consider is a simple model of 
two non-interacting non-relativistic free particles in one 
spatial dimension that has been "parameterized" , that 
is, Newtonian time is introduced as a canonical variable 
conjugate to the energy. The reader may question how 
relevant these simplified examples are to the issue of in- 
terest, namely theproblem of time in quantum gravity. 
To quote Kuchaf [l[ "The nature of the conditional prob- 
ability interpretation is so clear from these examples that 
it is hardly necessary to spell out how the formalism looks 
in quantum gravity". The reader will confirm this point 
of view while seeing how one gets the result for the prop- 
agator virtually without using any special features of the 
model in question. In particular, although the model 
does have a naturally defined time variable, we only use 
it to construct easily the evolving constants. The latter 
are known to exist in many examples (e.g. where 
there is no natural decomposition of the constraint into 
the "p° + H" form. 

The system has three configuration variables qo,qi,q2 
and the corresponding canonical momenta po,pi,P2- 
There is a constraint <\> = p +p\ /(2toi) +p 2 / {2ni2). The 
gauge invariant quantities, which have vanishing Poisson 
brackets with the constraint, are Q\ = q\ — piqo/mi and 
Q2 = 92 — P2qo/m2 and p\ and p2- These Dirac ob- 
servables represent the initial position and momenta of 
the particles. We then define evolving constants -^i(^) = 
Qi +Pit/m 1 and X 2 (t) = Q 2 +p2t/m 2 . We can check 
that they have vanishing Poisson bracket with the con- 
straint and that when t = qo then X\(t = qo) = Qi and 
X%(t = qo) = 172 • The quantization of the model is imme- 
diate [H . The states that are annihilated by the quantum 
version of the constraints are given by il>(pi,P2) times a 
prefactor 8(po + pf/(2mi) + p^/ {2rri2)) and the Hilbcrt 
space is that of square integrable functions ip(pi,p2), or 
equivalently in Fourier space by functions ip{q\ 1 q2). In 
this Hilbert space the evolving constants are well defined 
operators. Their common eigenstates are of the product 
form, 



1p Xl ,x 2 ;t = (j>l,P2\Xl,X 2 ;t) 



(3) 
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exp 



P1X1 +P2X2 - t 



( P\ 



p\ 



\2mi 2rri2 



with eigenvalues x±,x 2 for X\,X 2 , at some value 
of t. With these we can construct the pro- 
jectors that appear in the conditional probability, 
p xi(t) = J*** dzi J x oo dz 2 \zi,Z2;t)(z 1 ,Z2;t\ and simi- 
larly for P X2 (t) . The limits of integration Xi± corre- 
spond to Xi ± Axi/2 where Axi is introduced since 
if one is dealing with variables that have continuum 
spectrum, one cannot ask for "the probability that qi 



is a given value Xi, but rather within an interval of 
width Axi centered at such value. Let us consider a 
physical state given by a Gaussian for both variables 
centered at two distant phase-space points x\ and x\ 
and p\ and p\, e.g. p = \ip )(ip \, with V>o(?i,?2) = 
Ylj-i exp (— (qj — Xj) 2 /Aj + ip°qj) ■ We can then com- 
pute the quantity in the numerator of |T]), (notice that the 
denominator is just given by the numerator integrated in 
Q\ from —00 to 00, so for brevity we only show explicit 
calculations for the numerator), 



Num (P (ct^l^i, X2, x\) , po) = / dt'dt 



(4) 



X2+ 



x / dy2dz 2 (x' 2 ,t'\y2,t)tf 2tt (il>% it ) {z 2 ,t\x' 2 ,t') 



x / dyidzi(x' 1 ,t'\y 1 ,t)il;° ut (ip° ut )* {zi^x'^t') 



where "0 



= 



[wi,t\ip l ) for i — 1,2, and we have used 



the fact that the density matrix for this model is of di- 
rect product form, namely \ipo) — | V'o ' V'o ) ■ Usually one 
would like to consider systems with this property which 
implies that the system under study and the clock do 
not interact (we are choosing x\ as the clock variable). 
The interval Axi must be taken much larger than Ai the 
width of the Gaussian in the state in order for the mea- 
surement of the clock variable not to "destroy the state of 
the clock" . A measurement with more precision implies 
a faster loss of the (desired) classicality of the clock. In 
the case of X2 we assume we are studying a microscopic 
variable {m,2 <C rni), i.e. that is behaving quantum me- 
chanically, therefore we may and will assume Ax2 much 
smaller than the width of the Gaussian A2 to simplify the 
calculation of the integrals by substituting mean values. 
Carrying out the integrations explicitly, 

Num (P (xzlx'i, X2, x\) , po) ~ / dt'dt\(x 2 , t'\x2, t)\ 2 

Jo 

x\^ t \ 2 Axje Axi (xl-^-t-x^J \^ X ,A 2 (5) 

where is a rectangular function that is unity in the 
interval of width Ax\ around its argument and zero oth- 
erwise. We have assumed that evolution times are small 
such that the value of Ax\ does not change significantly. 
We have approximated the integrals in yi and z\ by in- 
tegrals from —00 to 00 since the Gaussian has a smaller 
support than the region of integration. The function 
arises since the approximation is good only if the peak of 
the Gaussian is within the integration region, otherwise 
the integral is close to zero. Putting together numerator 
and denominator we get, 



P ' (x 2 \x / 1 ,X2,X 1 ,p Q ) 



(6) 



lim / dt'\(x , 2,t'\x2,t(x 1 )}\ 2 V x '(t')Ax2 



4 



where t(x\) is the central value of t de- 
termined by the function and V x > (t') = 
Ti(P x[ (t')p ) / J^dtT^P^po) can be interpreted 
as the probability that the external (unobservable) time 
q° is t' when the variable we take as clock reads x[ . This 
would be controlled by the position of the peak and 
width of the Gaussian in the quantum state we chose. 
If instead of a Gaussian one had a Dirac delta, then we 
would recover the correct ordinary non-relativistic prop- 



agator, P(x' 2 \x' 1 ,x 2 ,x 1 ,p ) 



(x' 2 ,t'(x' 1 )\x2,t(xi)}A 2 , 



where t(x±) is determined by our choice of initial state 
to approximate the ordinary non-relativistic time corre- 
sponding to the position X\. The resulting expression is 
an approximation to the integral in x 2 of the ordinary 
propagator, therefore the factor £±x 2 . As is expected in 
relational treatments, one only obtains the traditional 
propagator at leading order. The use of real clocks leads 
to loss of quantum coherence, as is well known Q and 
therefore to corrections to the ordinary propagator. No- 
tice that up to now discussions of loss of coherence due 
to real clocks did not involve the presence of constraints, 
since they were framed for the gravitational case in the 
context of uniform and consistent discretizations Q, 
where constraints are eliminated. Here we confirm the 
presence of these effects in totally constrained systems. 

Let us sketch how the above proposal could be imple- 
mented in the case of general relativity. We consider 
the theory in vacuum coupled to a clock. We char- 
acterize the clock by its worldline X^{t) and T(t) its 
proper time. The action is the Einstein-Hilbert action 
for general relativity plus a term for the clock of the 

form S = -mj dr^J -X^^g^X^)) - f 2 where the 
dots mean total derivative with respect to the parame- 
ter r, and m is the mass of the clock. The equations 
of motion state that A M is a geodesic of the metric g^ 
and an equation stating that T is proportional to the 
proper time. As usual, we are assuming that the clock is 
a probe and therefore ignore back reaction. Classically 
this is certainly a good approximation. In this system we 
have only introduced a clock, not a complete coordinate 
system, one can ask only certain relational questions. For 
instance what is the value of a geometric scalar (e.g. the 
Kretschmann invariant) K at the space-time position of 
the clock when the clock measures a given value of time 
T. To complete our proposal one needs to find evolving 
constants S, functions of the metric and its first deriva- 
tives, parameterized with four real parameters x^ such 
that when they equal certain combinations of the metric 
and its derivatives the evolving constants reproduce the 
geometric quantity S we want to measure. The explicit 
construction of these quantities in general relativity can 
be onerous, but progress can be done by perturbative 
techniques, for example (see [l(|). One can then define 
the relational probabilities that the geometric quantity of 
interest take the value So when the clock measures time 
To, 

J d 4 xTr(^p)TT (P So (x)P To (x) P P To (x)) 



Defining a propagator needs more work, namely setting 
up a full coordinate system (i.e. introducing rulers in 
addition to clocks or considering a cloud of clocks as in 
[ill]). The calculational complexity would be important 
but modern loop quantum gravity techniques may allow 
a proper calculation. The expressions of the conditional 
probabilities in a situation like general relativity will not 
only include loss of coherence in time but also spatially, 
as has been analyzed in field theory in [lij ]. 



Summarizing, we have shown that one can formulate 
a completely relational picture of evolution in generally 
covariant systems framed entirely in the physical space of 
states and that yields the correct propagators in model 
systems and opens the possibility of assigning probabil- 
ities to histories and characterizing consistently the dy- 
namics of quantum general relativity. The resulting the- 
ory also predicts a fundamental mechanism of decoher- 
ence similarly as the one originally discussed in Q . 
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